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Determining g& using non-perturbatively 0(a) improved Wilson 
fermions 
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A completely non-perturbative estimate is given for qa using both quenched and unquenched 0(a) 
improved Wilson fermions. Particular attention is paid to the determination of the axial renormalisation 
constant, Za, using the Ward identity for the propagator. For the quenched case, we have results at three 
lattice spacings allowing a continuum extrapolation. 



1. INTRODUCTION 

The axial charge of the nucleon, Aq, is de- 
fined from the axial current by 

(N(p,<t)\A R \N(p,Z)) = WAq, 

with A R — (97^75 g),R, giving for the non- 
singlet charge 

g A = gf = Au- Ad. 

This determines the strength of /3-decay and 
also occurs in the Bjorken sum rule where 
the lowest moment of the difference between 
the proton and neutron polarised structure 
function g\ is proportional to the non-singlet 
charge. To find gA the programme is: 

1. Compute (N\uj fl "f^u-~d'y tIl ^d\N) on the 
lattice for many quark masses (3> 3) and 
many lattice spacings (3> 3). 

2. Determine non-perturbatively the renor- 
malisation constant Za- 

3. Chirally extrapolate m q (or m 2 ) — » 0. 

4. Continuum extrapolate a 2 — > (if using 
0{a) improved fermions) 

5. Do first for quenched (as much cheaper 
in CPU time) then repeat for un- 
quenched fermions. 

6. Compare with experiment, gA ~ 1.26. 



Our lattice simulation, point ([!]), is standard 
using ratios of nucleon three-point to two- 
point correlation functions and will not be dis- 
cussed further here, except to note that be- 
cause we are considering a non-singlet func- 
tion, the difficult to compute quark line dis- 
connected terms cancel. Here, in this talk, we 
shall mainly discuss point (|J). 

2. DETERMINING Z A 

On the lattice, for Wilson fermions, the chi- 
ral Ward Identity is only approximately true, 
due to discretisation effects. We have a fi- 
nite renormalisation constant ZA{g) to deter- 
mine by demanding that a Ward Identity is 
obeyed to 0(a n ) (with n = 2 for 0(a) im- 
proved fermions). The ALPHA collaboration 
has found Za for quenched fermions using 
the Schrodinger functional and PC AC; here, 
alternatively, we shall use the chiral Ward 
identity for the propagator. We shall first 
check our Za results for quenched fermions 
before then determining Za for unquenched 
fermions. The chiral Ward identity for the 
propagator, Sr, reads 

7 5 S^ 1 + S^js = 2m R A R , 

where ttir is the quark mass, and A^ the 
1PI pseudoscalar vertex, obtained from the 
Green's function Gp by amputating the exter- 
nal propagators. On the lattice we take this 
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equation to be correct to 0(a 2 ), 
TrS: 1 = Z A mTr l5 At 



where 0(a) improved operators are denoted 
with a star. The bare quark mass to* has also 
been replaced by the quark mass obtained^ 
from PC AC, rh. These quark masses may be 
computed in the standard way using the ap- 
propriate axial vector matrix elements. The 
propagator and vertex can also be found on 
the lattice, in the Landau gauge, A prob- 
lem with using the propagator is that we now 
need 0(a) improvement off-shell. To try to 
ameliorate this problem we note that as dom- 
inant off-shell effects come at short distances, 
we shall subtract a contact term (ie delta func- 
tion in position space), 0, from the propaga- 
tor and Green's function in addition to im- 
proving the appropriate operator. This proce- 
dure has been shown to work to first order in 
perturbation theory, j^j, and we shall use the 
formulation given there. Note however that 
the pseudoscalar Green's function (or IPI ver- 
tex) is particularly simple: there is effectively 
only one additional 0{a) improvement term. 
However, initial attempts to satisfy the above 
equation fail because of the presence of large 
0(ap 2 ) terms mainly due to the Wilson term 
in the propagator, and so here we choose to 
re-write the propagator and vertex Af in 
an 0(a) equivalent form, 0, as solutions of 



s- 1 



and 



1 — amh 



J» 5« A * on iJ* j 



[1 - amdp] Af = A p + iaX^S' 1 , Af }, 

(p M = (1/a) sinap^, p^ = (2/a) smap^/2) 
where by = -1 + 0(g 2 ) and Xy = 1 + 0(g 2 ) 
are the improvement coefficient and coefficient 
of the contact term for the propagator and 
dp = by - XpZ rn Z P + b A = -1 + 0(g 2 ). 
This gives S^ 1 = if + to* + 0(a?m 2 p 2 ) for 
the free field. The Wilson term has been sup- 
pressed, 0(ap 2 ) terms being multiplied by an 

A small additional factor 1 — am(bp — b A ) has been 
absorbed into a redefined A^. 



P=6.20 

Z A =0.804(10) 
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Figure 1. ThS" 1 plotted against mTry 5 A^ for j} = 
6.20 for 4 quark masses. Heavier masses lie on the 
RHS of the picture. Smaller momenta are on the RHS 
of each data set. 



additional small a 2 m 2 term. As in this case 
Tr 75 Af - 1/(1 + \ {am) 2 ) it is little effected 
from its true value of 1. 

We shall first determine Za for quenched 
fcrmions. We have generated data sets at /3 = 
6.0, 6.2, 6.4, each data set with 3 (or more) 
quark masses. To determine Za, we attempt 
to minimise R(p) = TiS~ x /fhTr^A^ = const 
for many momenta (a three parameter fit). A 
typical result is given in Fig. @. The gradient 
gives an estimate for Za- Deviations are seen 
for the higher momenta. This is more clearly 
seen in Fig. || which tells us how far we can 
go, a good fit being obtained for (ap) 2 < 2.5, 
say. Once the quark propagator has been ob- 
tained, we may then use it to find Z^ OM and 
m MOM ^ |j=j|_ ghgjj n0 ^ j-^ concerned with 



this here, g], but just concentrate on Za- In 
Fig. |§| we compare our result here with the 
previously known result, |lj. Reasonably good 
agreement is found. (In any case, as different 
methods have different 0(a 2 ) corrections com- 
plete agreement does not have to be found.) 

Buoyed up by this we now turn to un- 
quenched rif = 2 fermions where we have anal- 
ysed three UKQCD data sets, |J, (t3,K sea ) = 
(5.29,0.1340), (5.26,0.1345) and (5.20,0.1350). 
These are matched at an approximately con- 
stant rp value, corresponding to a ~ 0.105fm 
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Figure 2. ThS" 1 plotted with Z A mTry 5 A^ for f3 = 
6.20 for the 4 quark masses (the heavier quark masses 
are higher). 





Za\c a =-0.02 


Za c a = -0.08 


(5.29,0.1340) 
(5.26,0.1345) 
(5.20,0.1350) 


0.65(2) 
0.68(3) 
0.64(1) 


0.73(2) 
0.75(2) 
0.71(1) 



Table 1 

Preliminary values of Za obtained for unquenched 
fermions. 

(using rrj = 0.5fm). A typical result is shown 
in Fig. The preliminary values obtained are 
given in Table [l]. Unfortunately for construct- 
ing m, the axial current improvement coeffi- 
cient ca is unknown; we shall use here two 
values, the first ca ~ —0.02 corresponding to 
a tadpole improved result, while the second 
ca ~ —0.08 is the non-perturbative quenched 
value at j3 = 6.0. While this does not change 
the fit, ie the product Zatti, the 10% change 
in m directly produces a 10% change in Za- 

As a further check we have first computed 
/„■. After performing the chiral extrapolation, 
for the three data points given in Table 0, we 
find the results given in Fig. [|. Although a 
satisfactory result is obtained, one should not 
read too much into this at present, as apart 
from the uncertainty in Z a only three quark 
masses are used in the chiral extrapolation. 



Figure 3. Za versus g 2 (filled circles). Also shown is 
the ALPHA non-perturbative determination, ^j, to- 
gether with the first order perturbation theory result. 
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Figure 4. TrS t 1 plotted against mTry^A^ . 

3. RESULTS FOR g A 

We are now ready to present our results for 
gA- Again after the (linear) chiral extrapola- 
tion^] we plot our results as a function of a 2 . 
The current results are shown in Fig. |. All 
results, at present, lie too low in comparison 
with the experimental number. Even allow- 
ing for the uncertainty in the determination 



2 A typical chiral extrapolation for /3 = 6.0 quenched 
fermions is shown in Fig. 1 of [TL CI] . (Note that we take 
the unknown improvement coefficient b A = here. We 
have checked that in the chiral limit, as expected, this 
plays no role in the numerical value.) 
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Figure 5. rof n versus (a/ro) 2 , with the quenched re- 
sults (filled circles, the continuum extrapolated value 
being given by an open circle), together with the un- 
quenched results, (open circles). Note here that be- 
cause f^ are decreases as ca decreases while Za in- 
creases, changes in ca have little effect on the result. 
The experimental value is shown as a filled square. 
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Figure 6. gA versus (a/ro) 2 , with the quenched re- 
sults (filled circles, the continuum extrapolated value 
being given by an open circle), together with the un- 
quenched results, (open circles). The experimental 
value is shown as a filled square. 

of Za and hence 17,4 for unquenched fermions, 
it would still seem that the result is around 
the quenched result. (Note that also find 
a low value of gA-) Further details will be 
published elsewhere, 

4. CONCLUSIONS 

A completely non-perturbative evaluation 
of a nucleon matrix element needs 

• large physical boxes 

• good determination of Z 

• chiral extrapolation (delicate) 



• continuum extrapolation (delicate) 

How far have we got? We have concentrated 
on gA here and tried to develop a way of de- 
termining Za using the propagator and pseu- 
doscalar vertex. It is apparent though that we 
are only at the beginning, not only for the un- 
quenched results, but even for the quenched 
results many more quark mass and a values 
must be used to allow for reasonably reliable 
chiral and continuum extrapolations. We fi- 
nally note that the successful determination 
of gA is a real test of QCD as it is a rela- 
tively simple (nucleon) matrix element, well 
measured in experiment. 
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